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Abstract We study the combined effects of convection and radiative diffusion 
on the evolution of thin magnetic flux tubes in the solar interior. Radiative 
diffusion is the primary supplier of heat to convective motions in the lower 
convection zone, and it results in a heat input per unit volume of magnetic flux 
tubes that has been ignored by many previous thin flux tube studies. We use 
a thin flux tube model subject to convection taken from a rotating spherical 
shell of turbulent, solar-like convection as described by Weber, Fan, and Miesch 
(2011, Astrophys. J., 741 , 11 ; 2013, Solar Phys., 287 , 239), now taking into 
account the influence of radiative heating on 10 22 Mx flux tubes, corresponding 
to flux tubes of large active regions. Our simulations show that flux tubes of 
<60 kG subject to solar-like convective flows do not anchor in the overshoot 
region, but rather drift upward due to the increased buoyancy of the flux tube 
earlier in its evolution as a result of the inclusion of radiative diffusion. Flux 
tubes of magnetic Held strengths ranging from 15 kG to 100 kG have rise times 
of <0.2 years, and exhibit a Joy’s Law tilt-angle trend. Our results suggest that 
radiative heating is an effective mechanism by which flux tubes can escape from 
the stably stratified overshoot region, and that flux tubes do not necessarily need 
to be anchored in the overshoot region to produce emergence properties similar 
to those of active regions on the Sun. 
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1. Introduction 

Sunspots are observable manifestations of magnetism on the solar surface, thus 
providing a link to the deep-seated dynamo mechanism. For the Sun, the pre¬ 
vailing dynamo paradigm suggests that the toroidal magnetic field giving rise to 
the longitudinal bands of solar activity is amplified and stored in a thin shearing 
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region at the base of the convection zone called the tachocline ( e.g. I Gilman 


2000l : Charbonneaul . 2010 1. Buoyant loops of magnetic flux must then traverse 


the bulk of the turbulent convection zone, eventually emerging at the surface to 
form active regions. 

The thin flux tube (TFT) approximation has been developed by a number 


Defouwl. 1976|; Roberts and Webbl. 1978: 

Soruitjl. 

1981a b: Moreno-Insertd. 1986 

Ferriz-Mas. Schiissler. and Antoni. 19891: 

Chen el. 

19921: Ferriz-Mas and Schiissler 


tions based on the TFT approach have provided valuable insight on the evolution 
of ri s ing m agnetic loops th rough a quiescent solar convectio n zone (see review by 
Fanl . 2009 ). More recently, Weber. Fan, and Mieschl 1 201 lh (hereafter Article 1) 


and Weber. Fan, and Miesch ( 2013l l (hereafter Article 2) have extended the TFT 


model to include the effects of a rotating spherical shell of turbulent solar-like 
convective flows on the evolution of rising flux tubes. Most previous TFT studies 
assume that flux tubes evolve adiabati cally throughout the convection zone (ex¬ 
ceptions to t his are e.a. Fan and Fisherl (Il996h: Moreno-Insertis. Schiissler, and 


Glampedakis ( 2002T ); iRemoell ( 20031 )1. This is a valid assumption for the upper 
~2/3 of the solar convection zone. However, in the lower ssl/3 of the convection 
zone, closer to the radiative interior of the Sun, there is a significant non-zero 
divergence of radiative heat flux due to the deviation from radiative equilibrium. 
Radiative diffusion in the lower convection zone is the primary supplier of heat 
to large-scale global convective motions. It also results in a heat input per unit 
volume [ Q v } of magnetic flux tubes, which may have a substantial effect on their 
buoyancy, and hence their dynamic evol ution. 

Due to radiative heating. iFan and Fisher ( 1996 1 found that flux tubes of 10 21 — 
10 22 Mx rise through the convection zone in «two to four months, a comparable 
or shorter timescale compared to flux tube s allowed to evolve adiaba t ically , which 
have rise times of ftdwo to ten months ( Fan . Fisher, and Delucal . 119931 1. This 


shorter rise time is a result of radiative heating near the bottom of the convection 
zone, which increases the density deficit (i.e. buoyancy) of the flux tube at deeper 
depths, initiating more quickly their buoyant rise toward the surface. These 
emerging flux loops also show similar qualitative features to solar active regions 
such as tilt angles and morphological asymmetries, not significantly different 
from simulafioi^resultsjvhereJlm^ubesevolveadiabaticall^^ejf^Fan, Fisher, 
and Del uca. 19931: D’Silva and Choudhuri . 19931: Ca li gar i , Mor e no-Insertis, and 


Schiis sler, 19951 : Caligari. Schiissler. and Moreno-Insertil 1998l h Fan and Fisher 


( 19961 1 also find that flux tubes subject to radiative diffusion will rise quasi- 


statically (i.e. all forces closely balance) through the convective overshoot region. 
Using typical values for the overshoot region as computed by numerical solar- 
structure m odels, flux tubes take ssone year or less to emerge from the ove rshoot 
region (e.a. Ivan Ballegooiienl . 1982 : Fan and Fisheil . 19961 : Rempell . l2003ll . This 


short storage time compared to the Rill-year solar cycle may have significant 
implications for the solar dynamo mechanism. 

The purpose of this article is to study the combined effects of turbulent 
solar-like convection and heating due to radiative diffusion on the evolution 
of active-region-scale 10 22 Mx magnetic flux tubes. We begin in Section [2] with 
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a description of our simulation model and the modification to the TFT energy 
equation due to the addition of the non-adiabatic heating by radiative diffu¬ 
sion. In Section [3] we address how the addition of radiative diffusion to the 
energy equation alters the dynamic evolution of flux tubes subject to solar-like 
convective flows. Emergence properties of these flux tubes, specifically latitude 
of emergence and tilt-angle properties, will be compared in Sections 14.11 and 14.21 
respectively, to those of flux tubes that evolve adiabatically. We will mention the 
problem of flux storage in Section [5l and end with a summary and discussion in 
Section [G] 


2. Model Description 
2.1. The TFT+ASH Approach 


We model the evolution of magnetic flux tubes in the solar interior following the 
thin flux tube (TFT) approximation. The TFT equations are derived from the 
ideal MHD equations, operating under the assumption that the flux tube radius 
[a] is small compared to the total length [L], and other relevant length-scale 
variations such as the local pressure scale height and radius of curvature of the 
flux tube. All physical quantities of the flux tube are taken as averages over the 
cross-section, only varying spatially along the flux tube axis. The TFT equations 
that describe the evolution of each Lagrangian element of the one-dimensional 
flux tube are as follows (and similar for Articles 1 and 2): 

B 2 \ B 2 , 

8^) + ^ k 

( 1 ) 


du 

p Tt = 


-2p(f2 0 x v) - (p e - p)[g - L? 0 x (!?o x r)] +1 — 


-C A 


Pe\(v - «e)j.|(« - Ve)_L 
(7T$/B)V2 


d 

dt 




\pj p 

OS 


ldp = J_dp P T ^i 
p dt 7 p dt ad p dt ’ 

pRT 

V = -, 




( 2 ) 

( 3 ) 

( 4 ) 


B 2 

P+—=Pe, (5) 

where, r, v , B, p 1 p , T, which are functions of the time [f] and arc length [s] mea¬ 
sured along the tube, denote respectively the position, velocity, magnetic field 
strength, gas density, pressure, and temperature of a Lagrangian tube segment, 
l = dr/ds is the unit vector tangential to the flux tube, k = d 2 r/ds 2 is the tube’s 
curvature vector, subscript _L denotes the component perpendicular to the flux 
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tube, $ is the constant total flux of the tube, p e , p e , and p, which are functions 
of depth only, are respectively the pressure, density, and mean molecular weight 
of the surrounding external plasma, g is the gravitational acceleration that is 
a function of depth only, f?o is the angular velocity of the reference frame co¬ 
rotating with the Sun, with Ho = 2.7 x 10 ” 6 rad s -1 , Cd = 1 is the drag 
coefficient, 7 is the ratio of specific heats, S is the entropy per unit mass, V a d 
is the adiabatic temperature gradient, and v e (r,t) is a time dependent velocity 
field (relative to the rotating frame of reference) that impacts the dynamics of 
the thin flux tube through the drag force term (last term in Equation <[T|)). The 
term v e accounts for both the local convective flows and mean flows such as 
differential rotation. Heating due to radiative diffusion is introduced through 
the last term on the right-hand side of Equation ([3]), which we will describe in 
detail in Section [2721 

In the above equations, we do not introduce an explicit magnetic-diffusion 
or kinematic-viscosity term. The thin flux tube is untwisted (*.e. magnetic held 
lines do not twist about the flux tube axis), and is discretized with 800 uniformly 
spaced grid points along its arc length [s]. The numerical methods used to solve 
for the flux tube evolution as determined by the above s et of equations have 
been described in detail by Fan. Fisher, and Delucal (1993). 

For the stratification and thermodynamic properties of the external field-free 
plasma, namely p e , Pe, T e , g , and S = V c — V a d, where V c = dlnT e /dlnp e and 
V a d is the value of V e that one obtains by considering local adiabatic pertur¬ 
bations, we use Hiereference_solaimnode] i- (hereaftei^^odel i- S} - ^^ i Christensen- 
Dalsgaard et al. ( 1996l f. See Figure 1 in Article 1 for profiles of 7f. o„. p e , and <5 
used in our simulations. Here we define the base of the convection zone as r cz b = 
0.723R.Q (5.026xl0 10 cm), the radius in Model S where the plasma changes 
from sub-adiabatic (stably stratified, V e < V a d) to super-adiabatic (unstably 
stratified, V e > V a d)- Model S contains no region of convective penetration 
below r cz b- To mimic a convective overshoot region, we extend Model S below 
r cz b with a simple polytropic, sub-adiabatically stratified region that decreases 
from values of S = 0 to — 10 ~ 3 over a distance of ss 0.01 R 0 (ss 10° cm) downward 
from r cz b- In this region, it is assumed that overshooting convective motions are 
efficient enough to establish a sub-adiabatic stratification such that the temper¬ 
ature gradient remains closer to the adiabatic value rather than the radiative 
temperature gradient of the deep interior (see Section 15.11 for more details). 

What sets the thin flux tube simulations introduced in Articles 1 and 2, and 
continued here, apart from previous simulations is the inclusion into the drag 
force term of a time-dependent convective velocity field v e (r,t) relative to the 
rotating frame of reference. This three-dimensional global convection simulation 
is computed se paratel y usi ng the Anclastic Spherical Harmonic (ASH) code, 
as described by Miesch. Brun. and Toomre 1 2006 1. The ASH code solves the 
3D anelastic Navier-Stokes fluid equations using a pseudo-spectral method with 
both spherical harmonic and Chebyshev basis functions, explicitly resolving the 
largest scales of motion, while treating small turbulent eddies with sub-grid 
techniques. 

The computational domain for the ASH simulation extends from r = 0.69 Rq 
to r = 0.97 Rq (4.8x10 10 cm to 6.75xl0 10 cm), and the density contrast across 


SOLA: rad_heat_arxiv.tex; 30 March 2015; 0:27; p. 4 















Effects of Radiative Diffusion and Convection on Rising Flux Tubes 




Figure 1. (Left) Snapshot of ASH simulation convective radial velocities used in this article 
at a depth of 23 Mm in a Mollweide projection, with the dotted line representing the solar 
radius r = R©. Strong downflow lanes (purple) at the boundary of giant convective cells 
surround upflow regions (yellow). Also known as banana cells, the structures at low latitudes 
are rotationally aligned and propagate prograde. (Right) Angular velocity (with respect to 
the inertial frame), averaged over longitude and time (time interval 755 days). Color table 
saturates at the values indicated, with extrema ranging from 326 nHz to 468 nHz. 


the domain is «69, or 4.2 density scale heights. The Rayleigh number [i? a = 
gr 2 dAS/ (vnCp)\ is 5xl0 6 in the mid-convection zone, where d = r2 — r\ is the 
depth of the domain, Cp is the specific heat at constant pressure, and v and 
k are respectively the turbulent viscosity and thermal diffusivity of the ASH 
simulation. The Reynolds number [R e = v rnls d/v\ is of order 50, where v lms is 
the root-mean-square velocity relative to the rotating reference frame. For more 
specific details on the ASH simulation used here, see Articles 1 and 2. 

A typical giant-cell convection pattern of this ASH simulation at a depth of 23 
Mm below the solar surface is shown on the left in Figure [T] Broad upflow cells 
are surrounded by narrow downflow lanes, which can reach maximum downflow 
speeds of nearly 600 m s -1 at a mid-convection zone depth of about 86 Mm 
below the surface. These downflow lanes align preferentially with the rotation 
axis at low latitudes, reflecting the presence of so-called “banana cells”. The total 
angular velocity of the convection simulation n/2ir (with respect to the inertial 
frame) is solar-like, and decreases monotonically from rM 70 nHz at the Equator 
to ~330 nHz at the Poles, and exhibits nearly conical contours at mid-latitudes 
(see Figure [II ri ght) , as observed in the solar convection zone (e.a. Thompson 
et ai, 20031 ). The combined influence of density stratification and the Coriolis 
force induces anti-cyclonic vorticity in expanding upflows and cyclonic vorticity 
in contracting downflows. These effects yield a mean kinetic helicity density 
[iJk] that is negative in the northern hemisphere and positive in the southern 
hemisphere (see Article 1 for an image of the associated kinetic helicity). Such a 
helicity pattern is typigglfQiy^glglingj-gompressiblc convectioj^gx^Miesch and 
Toomre, 20091 ). 
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We emphasize that we use the one-dimensional Model S to describe the 
background stratification and thermodynamic properties of the computation 
domain, utilizing the ASH simulation only for the convective flow field. This 
choice was made to facilitate comparison between flux tubes allowed to evolve 
both with and without convection (see Section Differences in results will 
then be solely due to the presence of convective flows or the form of the TFT 
energy equation (see Section r2.2ll . In any case, the ASH background reference 
state is not significantly different from Model S, by at most «10% throughout 
the majority of the simulation domain. Time-varying temperature fluctuations 
in the ASH simulation from its prescribed reference state are small, on the order 
of S in the bulk of the convection zone. Furthermore, the ASH simulation is 
computed in the hydrodynamic regime separately from the TFT simulations. 
Any thermodynamic changes associated with the presence of the magnetic field 
in a truly three-dimensional magnetohydrodynamic ASH simulation would be 
on the order of 1//3, where /? = 8irp e /B 2 . 

Our simulations start with toroidal magnetic flux rings in mechanical equilib¬ 
rium located at a radial distance to the center of the Sun r = r 0 = 5.05 x 10 10 
cm, slightly above the base of the solar convection zone at r = r cz b = 5.026 x 10 10 
cm. At ro, the super-adiabaticity is S ss 3xl0~ 8 . To ensure the initial state of 
mechanical equilibrium, the flux tube is made neutrally buoyant by reducing the 
internal temperature compared to the external temperature and the magnetic 
tension force is balanced by the Coriolis force from an initially prograde toroidal 
flow in the tube. The external time-dependent convective velocity field described 
impacts the flux tube through its associated drag force (last term in Equation 
©)• 

We consider a range of flux tubes with initial magnetic field strengths [Ho] 
of 15, 30, 40, 50, 60, 80, and 100 kG, and initial latitudes [#o] ranging from 1° 
to 40° in either hemisphere. Considering that the root-mean-square (rms) of the 
convective downflows from the ASH simulation at the base of the convection 
zone is ~35 ms -1 , the equipartition magnetic field strength for our simulation 
is the B eq « 5 kG. In this case, we are investigating flux tubes on the order of 
3 — 20 B eq . The flux of the tube is constant at 10 22 Mx for this article, typical of 
large active regions with the strongest sunspots. As was also done in Articles 1 
and 2, each flux tube is initially perturbed with small undular motions consisting 
of a superposition of Fourier modes with azimuthal order ranging from to = 0 to 
to = 8 with random phase relations. However, such perturbations have negligible 
effects on the evolution of flux tubes subject to convection because perturbations 
provided by the convective velocity field are much stronger in amplitude. 

We perform seven groups of simulations sampling different time ranges of the 
ASH convective flow at each magnetic field strength and initial latitude in both 
the northern and southern hemispheres, for a total number of «2300 individual 
flux tube simulations. The flux tubes comprising one group are released at the 
base of the convection zone at the same starting time, although they do not 
interact with each other (i.e. are isolated) and are allowed to evolve until some 
portion of the flux tube reaches the top of the simulation domain. The flux tube 
release times for the groups are arbitrarily chosen, but are at least separated 
by the convective turnover timescale of the ASH convection simulation, which 
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is ~30 days. In this way, the flux tubes are able to sample significantly different 
portions of the convective velocity field. This procedure is the same as was done 
in Article 2. 

2.2. The Thin Flux Tube Energy Equation 

In the radiative interior of the Sun, energy is transported only by radiation. 
However, in the outer third of the Sun by radius, convection begins to take over 
much of the energy transport. Therefore, in the convection zone and overshoot 
region, the plasma is no longer in radiative equilibrium, resulting in a non-zero 
divergence of radiative heat flux. This acts to heat the flux tube in the lower 
convection zone, and cool it near the surface. Additionally, there is radiative 
diffusion across the flux tube due to a difference in temper ature between the 
interior of th e flux tub e and the surrounding pl asma fscc. e.o.. |Parkerl.ll979l: van 
Ballegooiien. Il982l h As a result, especially in the lower convection zone, it does 
not suffice to assume that the flux tube evolves adiabatically. The incorporation 
of radiati ve heating to the TFT equations has been described in detail bv Fan 
and Fisher" (fl996t l. Here we briefly review their methods for clarity in this article. 
The rate of heat in put per unit volume [dQ v /dt] of the flux tube plasma is 


([Fan and Fisher, 1996j): 


dS dQ v _ , 

pT Ht = nr = v ' (Kvn 


( 6 ) 


where S is the entropy per unit mass, the right-most expression is a divergence 
of a radiative diffusive heat flux [kVT], and n is the coefficient of radiative con¬ 
ductivity. Note that in the adiabatic regime dS'/dt = 0. This is the assumption 
made in Art icles 1 and 2. 

Following iFan and Fisher ( 1996t h Equation © can be reduced to two domi¬ 
nant terms: 

d Q 


df 


V ■ (AC e VT e ) + kV]_(5T, 


(7) 


where K e and k are the coefficient of radiative conductivity of the background 
field-free plasma and the flux tube plasma, respectively. The term ST = T — 
T e represents a small difference between the external and internal temperature 
introduced due to the presence of a magnetic field in the flux tube, and Vj_ is 
the gradient vector in the direction perpendicular to the flux tube axis (in the 
plane of the flux tube cross-section). The first term on the right-hand side of 
Equation o, which we will subsequently refer to as (dQ/dt)\ : is dependent only 
on the external plasma thermodynamic quantities n e and T e , which vary only 
as a function of radial distance from Sun center. This term is zero if the fluid is 
in radiative equilibrium, which is the case in the radiative interior of the Sun. 
The second term, subsequently referred to as (dQ/dt) 2 , represents a radiative 
diffusion across the flux tube due to the temperature difference between the flux 
tube and the ex ternal plasma. The two heating terms are computed as follows 
( Fan and Fisher . 1996i ): 


dQ 

dt 


= -F tc 


dr VV rad 


( 8 ) 
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(9) 


where the magnitude of the total energy flux is given by F to t = L/(47rr 2 ), and L 
is the total luminosity of the Sun, V c = dlnT e /dlnp e is the local temperature 
gradient of the external plasma environment, and V ra d = (dlnT e /dlnp e ) ra( j 
is the temperature gradient required for the energy to be transported only by 
radiation, on ss 2.4048 is the first zero of the Bessel function Jq(x), a is the flux 
tube radius, and ST = T — T e is the mean temperature difference across the flux 
tube. 

The magnitude of (dQ/df)i ~ n c X7 c T e /H pl where H p = p e /gp e is the pressure 
scale height and we have used d/dr = 1/H„ and /H„ ( Fan and 


Fisher, 19961 b For a neutrally buoyant flux tube at the base of the convection zone 
p = p e . Such that the condition of pressure balance (Equation (|5]l) is satisfied, 
the internal temperature of the flux tube must be less than the surrounding 
environment. This temperature deficit (5T/T e ) ss 1//3 = B 2 /8np e , where /? 1 
in the solar interior. Then, comparing the magnitudes of (dQ/d t)i and (dQ / dt) 2 
to de termine their relative importance to the energy equation (IFan and Fishcrl . 
19961 ): 


Ve 0^7; 

a 2 H 2 ST 


8 V e Pe $ 
a? HlB 3 ' 


( 10 ) 


Throughout the convection zone V c ~ 0.4. Taken from Model S, typical values 
of the pressure scale height and pressure at r 0 = 5.05 x 10 10 cm near the base 
of the convection zone are H p « 5.5xlO 9 cm and p e ~ 4.7xl0 13 g s -2 cm -1 . 
Therefore, the magnitude of Equation (1101) is «8.6xl0~ 7 ®/B 3 . For $ = 10 22 
Mx, the ratios of |(dQ/df)i| to |(dQ/dt) 2 | following this relationship are given 
in Table [1] As Equation (fTUl) is proportional to $, the values quoted in Table [T] 
will decrease by an order of magnitude for $ = 10 21 Mx, and by two orders of 
magnitude for <f> = 10 2 ° Mx. 

The second heating term [(dQ/df) 2 ] acts to reduce the temperature difference 
between the flux tube and the external plasma environment, bringing it closer 
to a state of thermal equilibrium, enhancing the density deficit of the flux tube 
in the process. This term is inversely proportional to the flux tube radius, there¬ 
fore it will increase in magnitude for thinner tubes (be. as the magnetic field 
strength increases and magnetic flux decreases). However, the first heating term 
[(d<3/di)i] is dependent only on properties of the background plasma as a func¬ 
tion of distance [r] from Sun center. Therefore each flux tube will experience the 
same heating from the term (dQ/dt)i at the same distance [r]. We have chosen 
in this article to focus on how radiative heating in conjunction with convection 
influences the dynamic properties of 10 22 Mx flux tubes, corresponding to large 
solar active regions. Following the values reported in Table [Q it is a reasonable 
approximation to neglect the term (dQ/di) 2 in the TFT energy eq uation for 
flux tubes of 10 22 Mx in the 15-100 kG range, as was done by I Fan and Fisher 
( 19961) . Future TFT studies incorporating the effects of radiative heating on flux 
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Table 1. Ratio of |(dQ/dt)i | to \(dQ/dt) 2 1 following Equation (1101) for flux tubes 
of 10 22 Mx. The term |(dQ/dt)i| is much larger than |(dQ/d£) 2 | for tubes of 10 22 
Mx. As |(dQ/d£)i|/|(dQ/di) 2 | is directly proportional to <l>, the quoted values 
will decrease by an order of magnitude for <f> = 10 21 Mx, and by two orders of 
magnitude for <f> = 10 20 Mx. Values have been rounded to two significant digits. 


15 kG 

30 kG 

40 kG 

50 kG 

60 kG 

80 kG 

100 kG 

|dQi/dQ 2 | 2500 

320 

130 

69 

40 

17 

8.6 


tubes of < 10 22 Mx will require the inclusion of (dQ/df) 2 , especially at larger 
magnetic field strengths. Combining the general TFT energy equation (Equation 
(HI wit h Equat i on ( 151 ) . the new energy equation for the TFT model we use is 
( Fan and Fishcrl . Il996f ): 


1 dp 
p dt 


1 dp 
7 p df 


V ad - 


F, r 


Y—y 

p dr VV rad y 


( 11 ) 


where the last term was assumed to be zero for the adiabatically evolving TFT 
simulations of Articles 1 and 2. Values for the quantities V ad , V e , V rad , and F tot 
are provided by Model S. 

This added heating term (last term of Equation (fTTIl l affects the dynamic 
evolution of the flux tube by changing its buoyancy evolution [dAp/df]. To 
illustrate the im portance of including the radiative-heating term in the TFT 
energy equation. iFan and Fisher! ( 1996 1 examined the growth of buoyancy caused 
by radiative heating as compared to that due to the adiabatic expansion of the 
flux tube plasma rising through a super-adiabatically stratified environment: 


dt J rad Pe \ dt 


( 12 ) 


dAp\ v 


dt ) 


ad 




4 A 


(13) 


where <5 = V e — V ad an d v T is the radia l veloc ity of the flux tube apex. It 
is found (see Figure 4 in (Fan and Fisherl (1996!)) that in the lower one-third 
of the convection zone, dAp/df is dominated by contributions from radiative 
heating. This term is only dependent on the properties of the background plasma 
environment, and is therefore the same for every flux tube. In the remaining two- 
thirds of the convection zone above «0.80 Rq, the evolution of the flux tube can 
be described as essentially adiabatic, as the heating rate [(dQ/dt)i] decreases 
with height and S and v T increase with height. As will be shown later in Section 
[3l the inclusion of the radiative heating term to the TFT energy equation has 
significant implications for the dynamic evolution and rise times of the flux tube. 

In Sections [3] and [H we will compare the results of the simulations calculated 
in Articles 1 and 2, where we assume that the flux tube evolves adiabatically, 
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Figure 2. Polar (top) and equatorial (bottom) view of flux tubes once some portion of the tube 
has reached the simulation upper boundary. These 10 22 Mx, do = 8° flux tubes are allowed to 
evolve in the absence of convection, but with the addition of heating due to radiative diffusion. 
The orange sphere has a radius of 4.9xl0 10 cm. The addition of radiative heating forces flux 
tubes to drift upward as a whole and away from the overshoot region. Only 100 kG flux tubes 
of 8o < 8° have footpoints that anchor in the overshoot region. 


with the new simulations as described here in Section [2j For simplicity, we will 
often refer to the simulations of Articles 1 and 2 as case AB (adiabatic). We will 
refer to the non-adiabatic simulations with the addition of radiative diffusion to 
the TFT energy equation as case RD. 


3. Flux Tube Dynamics 

3.1. Flux Tube Morphology 

When flux tubes from our simulations evolve adiabatically in the absence of 
convection, rising buoyant loops develop solely as a result of the non-linear 
growth of the magnetic buoyancy instability, as discussed in Article 1. The 
troughs of these rising loops penetrate into the overshoot region where they 
remain anchored for the duration of the flux emergence process. However, when 
heating due to radiative diffusion is considered, only 100 kG flux tubes of 9q < 8° 
are capable of anchoring in the overshoot region (see Figure [5]). 

The term (dQ/df)i acts to heat the flux tube uniformly (i.e. uniformly for 
flux tube portions at the same distance [r]), especially near the base of the 
convection zone where the divergence of radiative heat flux from the external 
plasma environment is the greatest. This uniform heating increases the flux 
tube’s density deficit early in its rise, thereby increasing the buoyancy of the 
flux tube in the lower convection zone. Only 100 kG flux tubes of 9 < 8° 
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Figure 3. Polar (top) and equatorial (bottom) view of flux tubes once some portion of the 
tube has reached the simulation upper boundary. These 10 22 Mx, Qq — 8° flux tubes are the 
same as in Figure [21 except the flux tube is subjected to the external convective flow. In all 
cases, the image has been rotated such that the flux tube apex is on the right, and at the three 
o'clock position if looking down from the north solar Pole. These flux tubes are referred to 
here as case RD flux tubes, indicating that radiative diffusion and convection are included in 
the simulations. Convection prevents the severe poleward slippage of the flux tube as depicted 
in Figure |2] and weaker magnetic-field-strength flux tubes are more susceptible to deformation 
by convection. 


develop undular magnetic buoyancy instabilities that grow fast enough to allow 
the troughs of at least one of their rising loops to penetrate into the overshoot 
region. Magnetic buoyancy instabilities develop more slowly in weak magnetic- 
field-strength flux tubes. Flux tubes of <60 kG develop a uniform buoyancy 
quicker than the growth of the magnetic buoyancy instabilities, resulting in a 
nearly axisymmetric rise. This is shown especially well for a 15 and 40 kG flux 
tube in Figure [2] 

When solar-like convective flows are included in the thin flux tube simulations 
(denoted as case RD flux tubes here because of the inclusion of radiative diffusion 
to the model), we find that 15 kG flux tubes are affected most strongly by convec¬ 
tion, whereas 100 kG flux tubes only have moderate deformations resulting from 
the strongest convective downdrafts (see Figure [3] and Figure 0 ]). This trend is 
the same as observed for case AB [i.e. flux tube simulations with convection, but 
only adiabatic evolution) in Article 1. The addition of convection also keeps flux 
tubes from suffering the severe poleward slippage shown in Figure [2] especially 
at weaker magnetic-field strengths. 

Figure 0 ] clearly shows that when convection is included, case RD flux tubes 
of <60 kG do not anchor in the overshoot region. Rather, due to heating from 
radiative diffusion, these flux tubes develop a density deficit very early in their 
evolution and begin to float away from the base of the convection zone before 
magnetic buoyancy instabilities of sufficient amplitude can set in to anchor the 
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Figure 4. Flux tube radial distance from Sun center r (black line), plotted with the external 
radial velocity experienced by the flux tube at the height r of the flux tube segment (heavy 
red line), both as functions of the azimuthal angle [</>]. Snapshots are for case RD flux tubes 
with <f> = 10 22 Mx, dg = 15°, for Bo decreasing from top to bottom. All flux tubes shown are 
initialized at the same time, and therefore experience the same flow field initially. The dashed 
line represents the base of the convection zone, below which is the stably stratified overshoot 
region. These plots show the evolution of the flux tube once the apex has reached a height 
of (left) 0.80 Rq, (middle) 0.88 Rq, and (right) 0.97 Rq, indicating that <60 kG flux tubes 
subject to radiative heating do not anchor in the overshoot region. 


flux tube in the overshoot region. They reach the middle of the convection zone 
where they are continually buffeted by convection until a buoyant loop reaches 
the simulation upper boundary. All flux tubes of >80 kG do anchor when con¬ 
vection is included, and their footpoints continue to drift out of the overshoot 
region as the tube evolves due to the uniform heating supplied by (dQ/dt)i. 
This anchoring is facilitated by perturbations to the flux tube provided by the 
convective flow field, initiating magnetic buoyancy instabilities of large enough 
amplitude to anchor the flux tube footpoints before the tube as a whole can 
drift away from the lower convection zone. These convective perturbations are 
larger in amplitude than the superposition of Fourier modes applied to the flux 
tubes to promote undular magnetic buoyancy instabilities when convection is 
not present. 

It is worthwhile to briefly mention the effects of mean flows from our ASH 
convection simulation on flux tube evolution. As mentioned in Article 1, differen¬ 
tial rotation present in the convection simulation tends to boost the azimuthal 
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Figure 5. The right-hand side of Equation (1141 ) (colored symbols) plotted for flux tubes 
evolving in a convective velocity field assuming adiabatic evolution (case AB, left), and with the 
addition of radiative heating (case RD, right). The right-hand side of Equation <11411 decreases 
with decreasing magnetic field strength, therefore the top symbol curve represents 100 kG flux 
tubes (diamond), followed by 60 kG (plus sign) and 40 kG (square), with the bottom symbol 
curve representing 15 kG flux tubes (cross). Also plotted are representative convective velocity 
field speeds in the radial direction from the ASH simulation (black lines). Flux tubes of <60 
kG are significantly more susceptible to convective influences when heating due to radiative 
diffusion is included. 


velocity of the flux tube mass elements in the prograde direction. We do not 
expect this to have a significant effect on the flux tube properties discussed in 
this article. The kinetic energy of the mean meridional flow present in our ASH 
simulation is approximately two orders of magnitude smaller than the kinetic 
energy of convection and differential rotation. For the range of flux tubes that 
we consider in this study, it is unlikely that meridional circulation contributes 
much to flux tube evolution. 


3.2. Convection vs. Magnetic Buoyancy 


As in Article 1, we compare the magnitude of the drag force to the magnetic 
buoyancy force acting on the flux tube to understand their relative importance 
on flux tube evolution for case RD simulations. Following Equations (6) and (7) 
from Article 1, for the drag force to dominate the buoyancy force in the radial 
direction: 


^cr 



(14) 


i.e. the convective flow speed [« cr ] needs to be greater than the Alfven speed 
[va. = B/(ATtpe) 1 / 2 ] multiplied by the term ( a/H p j 1 / 2 . 

In Figure [5J we have plotted the right-hand side of Equation (1141) for flux 
tubes of $ = 10 22 Mx and do = 15° for four different initial magnetic field 
strengths. For comparison, the left panel shows the right-hand side of Equation 
(ED for case AB, and the right panel for case RD. Each flux tube is subjected 
to the same flow held, and representative radial convective-flow speeds are also 
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Figure 6. (Left) AT = T e —T and (Right) magnetic pressure at the apex of a flux tube as a 
function of height, where = 10 22 Mx, Bo = 60 kG, and Oo = 15°. These quantities are shown 
for a case AB flux tube (plus symbols), and a case RD flux tube (diamond symbols), both 
rising through a turbulent solar-like convection zone. Heating supplied by radiative diffusion 
increases the internal temperature of case RD flux tubes as compared to case AB flux tubes, 
resulting in a comparative increase in internal gas pressure which contributes to a decrease of 
magnetic pressure. 


plotted. Flux tubes of B 0 = 100 kG in both the left and right panels of Figure 
[5] are only affected by the strongest downflows, and the plots of v a (a/Hp) 1 ' 2 at 
the apex of the loop are very similar. However, at weaker initial magnetic field 
strengths, the quantity v a (a/Hp) 1 / 2 is decreased in the upper convection zone for 
case RD flux tubes as compared to case AB flux tubes. Figure [5] shows that case 
RD flux tubes of <60 kG are affected more strongly by convection than case AB 
flux tubes of the same initial magnetic field strength. This occurs because the 
field strength of the flux tube apex becomes weaker for the case with radiative 
heating. 

Due to radiative diffusion in the lower convection zone, flux tubes experience 
an increase in their internal temperature. The quantity AT = T e — T at the apex 
of a case AB and case RD flux tube, both subject to the same convective flows, 
is shown in the left-hand side of Figure [6] In the lower portion of the convection 
zone, heating due to radiative diffusion increases the internal temperature of case 
RD flux tubes as compared to case AB flux tubes. This results in a comparative 
increase in the internal gas pressure of case RD flux tubes, contributing to a 
decrease in the magnetic field strength so that the condition of pressure balance 
is fulfilled (Equation ( 0 ). The magnetic pressure at the apex of a case AB and 
case RD flux tube is also shown in the right-hand side of Figure [6] Indeed, the 
magnetic pressure, and therefore the magnetic field strength at the apex of the 
case RD flux tube is less than it is for the case AB flux tube throughout the 
bulk of the convection zone. This reduced magnetic field at the flux tube apex 
implies that the flux tube will be advected more strongly by convection. 

3.3. Rise Times 

Adiabatically evolving flux tubes spend the majority of their rise times at the 
base of the convection zone as the magnetic buoyancy instability grows. In our 


SOLA: rad_heat_arxiv.tex; 30 March 2015; 0:27; p. 14 








Effects of Radiative Diffusion and Convection on Rising Flux Tubes 


0.8 


0.6 


0.4 


0.2 


0.0 


X Adiabatic 
\> Radiativb 
Diffusion - 


X 


0 


<> 


X 

❖ 


o 


20 40 60 80 100 

Initial Magnetic Field B (kG) 

Figure 7. Average rise times for <3? = 10 22 Mx case AB (crosses), and case RD (diamonds) 
flux tubes. Radiative heating reduces the average rise time of the flux tube in most cases. 


simulations, convective flows provide perturbations to the flux tube that initiate 
the growth of these instabilities. However, when radiative heating is included, 
the flux tube receives a kick start in its rise toward the surface, as the flux tube 
is heated uniformly by the non-zero divergence of radiative heat flux in the lower 
convection zone, thereby increasing its density deficit [A p = p e — p\ earlier on, 
speeding up its rise through the solar convection zone. 

The inclusion of radiative diffusion in the TFT energy equation reduces the 
flux tube rise times to <0.2 years for all magnetic field strengths, as compared 
to <0.7 years for flux tubes evolving adiabatically (see Figure[7]). The reduction 
is minimal for 80 and 100 kG flux tubes. In this large-magnetic-field-strength 
regime, the average rise times are nearly the same because of the strong magnetic 
buoyancy and tension of the flux tubes, and because their undular magnetic 
buoyancy instabilities of order m = 1 to m = 3 grow quickly, facilitating 
anchoring in the overshoot region. 

There is a large reduction of rise times for case RD flux tubes in the Bq = 
15 — 60 kG regime as compared to case AB flux tubes: a difference of «0.05 — 0.5 
years. This reduction is a result of radiative heating contributing a significant 
increase to the buoyancy of the flux tube early in its rise. Radiative diffusion 
also increases the internal pressure of the tube early on, forcing a reduction in 
its magnetic field strength, therefore it is more susceptible to convection. 

Case AB flux tubes also show a large spread in their average rise times, 
from «0.15 — 0.7 years. This is a result of the complex interaction between the 
buoyancy and drag forces on the flux tube for various magnetic field strengths 
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(see Articles 1 and 2). Case RD flux tubes only exhibit a spread in rise times 
of ~0.1 — 0.2 years. The average rise times are brought closer together for all 
magnetic-field-strength flux tubes because they all experience the same heating 
rate at the same height in the deep convection zone. 


4. Emergence Properties 
4.1. Latitude of Emergence 

In the absence of convection, a flux tube will emerge radially with 9 em ~ 
(emergence latitude ~ initial latitude) if the outward component (away from 
the rotation axis) of the buoyancy force domina tes the inward component (to - 
ward the rotation axis) of the Coriolis force (e.a. IChoudhuri and Gihnanl . ll987 ). 
Such a scenario occurs for flux tubes of large magnetic field strength, i.e. 100 
kG (see Articles 1 and 2). Conversely, a flux tube rises more parallel to the 
rotation axis as the outward component of the buoyancy force is balanced by 
the inward component of the Coriolis force. This happens as the initial magnetic 
field strength of the flux tube is reduced, i.e. B 0 < 60 kG, leading to a poleward 
deflection of the flux tube away from its initial latitude [$o] that increases as Bq 
decreases (see also Articles 1 and 2). 

When convection is introduced, case AB (i.e. adiabatic.ally evolving) flux 
tubes of large B 0 , and therefore larger buoyancy, still rise mostly radially (left 
panel of Figure 0 ). However, especially at weak field strengths of <30 kG, the 
buoyancy force is reduced compared to the Coriolis force acting on the flux tube, 
forcing the apex of the flux tube to deflect poleward (i.e. more parallel to the 
rotation axis). This effect is responsible for the moderate latitudinal deflection 
of 15 kG flux tubes, especially at high latitudes (see Figured lower left panel). 
As shown in Articles 1 and 2, the addition of solar-like convective flows prevents 
weak magnetic-field-strength flux tubes from deflecting poleward as severely 
as when convection is not included. Convection also introduces a scatter in 
the emergence latitude values, which increases as Bq decreases. This happens 
because flux tubes become more susceptible to advection by convective flows as 
both the magnetic buoyancy and tension forces are reduced. 

Case RD (i.e. radiative diffusion included) flux tubes of 80 — 100 kG emerge 
with latitudes similar to case AB flux tubes of the same field strength (right panel 
of Figure [8]). This is not surprising, as flux tubes in this field-strength regime 
do anchor in the overshoot region. However, as the magnetic field strength of 
these flux tubes decreases, we begin to notice that case RD flux tubes tend to 
emerge at latitudes larger than case AB flux tubes, especially in mid-latitudes. 
Due to the conservation of angular momentum, as case RD flux tubes rise as a 
whole away from the convection zone base, a retrograde flow of plasma inside the 
flux tube is enhanced, increasing the Coriolis force acting on the now rising flux 
ring, deflecting the case RD flux tubes more poleward. At low latitudes for case 
RD flux tubes of <60 kG, convective flows are able to keep the flux tube from 
slipping poleward, and there are no latitudinal zones void of flux emergence, as 
shown for case AB flux tubes in Article 1. 
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Figure 8. Initial latitude versus emergence latitude of the apex of 10 22 Mx flux tubes. The 
left column shows the latitudinal deflection for case AB flux tubes, and the right column 
for case RD flux tubes (labeled as non-adiabatic here). Both axes are in units of degrees. 
A dashed line indicates where the emergence latitude equals the initial latitude. With the 
addition of radiative heating, middle-to-high latitude flux tubes of moderate to weak magnetic 
field strength deflect poleward. 
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4.2. Tilt Angles 


Solar active regions tend to emerge with their leading polarity (in the direction 
of solar rotation) closer to the Equator than the following, such that a line drawn 
between the center o f the two bipolar regions will be tilted with respect to the 
East-West direction ( Hale et al. . 1919 ). Here, the tilt angle is computed as the 
angle between the tangent vector at the apex of the emerging flux tube (once it 
has reached the top of the simulation domain), and the local East-West direction. 
We define a positive sign of tilt as a clockwise/counter-clockwise rotation of the 
tangent vector away from the East-West direction in the Northern/Southern 
hemisphere, consistent with the direction of the observed mean tilt of active 
regions. If the magn itude of the tilt angle exceeds 90°, then the active region 

19191) . arity in 


violates Hale’s Law (Hale et al. 


the direction of solar rotation for that particular hemisphere. 

To align our statistical results with those obtained from observations, unless 
otherwise stated, tilt angles derived from our flux tube simulations that do not 
fall in the range [-90°, 90°] are shifted to be brought back into this interval, 
thereby losing information about anti-Hale tilt angles. This approach is usually 
what is done for tilt-an gle st atistical studies of active regions on the Sun, starting 
withlHale et al. (fl919l) (e.q. Wang and Sheelev . 1989 ; Fisher. Fan, and Howardl 


1995 : Howardl. 199fil: iDasi-Esmiig et al. 201fll: Stenflo and Kosovichev. 20121: Li 


and Ulrich. 120121 : McClintock and Nortonl . 120131 ). If the leading polarity of the 

bipolar region in the direction of solar rotation, regardless of its sign, is closer 
to/farther from the Equator than the following, it will then be identified as a 
positive/negative tilt in the range [0°, 90°]/[-90°, 0°]. This approach is different 
from what was done in Articles 1 and 2, where tilt angles were in the range of 
[-180°, 180°]. In retrospect, we feel that shifting the tilt angles to fall in the [-90°, 
90°] range is probably a better diagnostic for comparison with observations. 

We perform numerous diagnostics on the tilt angles of our simulated flux 
tubes in an attempt to better constrain the magnetic field strength at which the 
solar dynamo may be operating, as was the topic of Article 2. Also, we wish 
to identify the effect radiative diffusion has on tilt angle trends, comparing the 
results of the 10 22 Mx case RD flux tube simulations calculated for this article 
with the case AB flux tubes from Articles 1 and 2. 


4-2.1. The Joy’s Law Trend 


Known as Joy’s Law, the average tilting behavior of emerging-flux r egions tends 


to increase in magnitude as the latitude of emergence increases ( e.g. IHalc et al 


Il919l ). Many authors have recovered Joy’s Law from observations, however there 
is generally no agreed upon common me t hod to obtain an empirical Joy ’s Law 


equation (e.g. Fisher. F an, and Howardl. Il995l: Dasi-Espuig et al . 2010 

uT 


flo and Kosovichev, 20121 : Li and IJlric a 120121 : iMcClintock and Nortonl . I 2 O Kil l. 


Sten- 


Within this section, we will employ two such methods in order to compare our 
simulation results to those of the Joy’s Law trend recovered from observations. 

For our first method (Method 1) of obtaining an empirical equation for the 
Joy’s Law trend from our simulations, we assume that the tilt angle increases 
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Table 2. Joy’s Law best-fit line slopes and rj \\ ti values for cases AB (Columns 2, 4, 6) 
and RD (Columns 3, 5, 7). (Column 1) magnetic field strengths. (Column 2, 3) Slopes 
m A (unitless) of the best-fit line following Method 1, and (Column 4, 5) slopes me 
(units of degrees) of the best-fit line following Method 2. The uncertainties reported 
in relation to the slopes are uncertainties in the determination of the fit parameters 
7 tia and mj, which takes into account both propagation of error and the variance of 
the data around the best-fit line. Slopes of the best-fit lines following Methods 1 and 2 
show a similar trend, reaching a maximum at 40 kG for case AB flux tubes, and 60 kG 
case RD flux tubes. (Column 6, 7) Standard deviation erg t (units of degrees) of the tilt 
angle about the best-fit Joy’s Law equation (Method 1). For magnetic field strengths 
of 40 — 60 kG, the values of ag t are substantially larger for case RD as compared to 
AB. 


B [kG] 

m A _ 

AB 

m A - 

RD 

m B 

i—AB 

m B 

i—RD 

CAB 

CRD 

100 

0.25 

± 

0.02 

0.32 

± 

0.02 

15.5° 

± 

i.i° 

19.8° 

± 

1.0° 

8.2° 

7.6° 

80 

0.29 

± 

0.02 

0.40 

± 

0.02 

17.7° 

± 

1.2° 

25.7° 

± 

1.2° 

9.1° 

10° 

60 

0.35 

± 

0.02 

0.44 

± 

0.04 

21.6° 

± 

1.5° 

27.1° 

± 

2.4° 

13° 

24° 

50 

0.38 

± 

0.03 

0.42 

± 

0.04 

23.8° 

± 

2.1° 

26.2° 

± 

2.6° 

18° 

26° 

40 

0.40 

± 

0.04 

0.27 

± 

0.05 

25.2° 

± 

2.5° 

16.7° 

± 

3.4° 

22° 

35° 

30 

0.31 

± 

0.06 

0.32 

± 

0.05 

19.6° 

± 

3.5° 

20.0° 

± 

3.3° 

33° 

34° 

15 

0.25 

± 

0.07 

0.31 

± 

0.07 

15.4° 

± 

4.2° 

19.1° 

± 

4.3° 

39° 

43° 


linearly with increasing emergence latitude. We perform a linear least-squares fit 
of the tilt angle as a function of emergence latitude following a = m A 9, where a 
and 8 represent the tilt angle and emergence latitude respectively, both in units 
of degrees, of our simulated flux tubes once the apex has reached the simulation 
upper boundary, and toa is the slope (unitless) of the best-fit line. The slopes 
toa of these best-fit lines along with their uncertainties are reported in columns 2 
and 3 of Table [2] Note that the values of toa for case AB flux tubes are different 
from those reported in Article 2 because of our choice here to restrict the tilt 
angles to [-90°, 90°]. The uncertainties on the slopes increase with decreasing 
magnetic field, as the flux tube is more susceptible to deformation by convective 
flows, and the tilt angles exhibit more of a scatter about the best-fit line. From 


white- light sunspot-group data spanning Solar Cycles 15 - 21. lDasi-Espuig et al. 


( 20101 1 find an empirical Joy’s Law equation of slope toa = 0.26±0.05 for Mount 
Wilson sunspot data, and toa = 0.28 ± 0.06 for Ko daikanal data. T he values 
we report in columns 2 and 3 of Table [2] agree with Dasi-Espuig et IzjO (120101) 
within the reported uncertainties for both cases AB and RD in the low (15-30 
kG) and high (100 kG) magnetic-field-strength regimes. 

A Joy’s Law fit can also be performed using the equation (Method 2) a = 
TOesind, which is a good choice assuming the origin of the tilt angle is related 
to the Coriolis force, as this force varies with latitude as sin(d). Here, the tilt 
angle a is again in units of degrees and tob also has units of degrees. The 
va lues we obtain from our s imula tions are shown in columns 4 and 5 of Table 
[5J IStenflo and Kosovichevl ( 20121) perform such a fit using 15 years of MDI 
full-disk magnetograms, finding a slope of tob = 32.1° ± 0.7 °. However, us ing 
Mount Wilson sunspot-group data from 1917 1985, Fisher. Fan, and Howard 
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(119951 ) find a best-fit equation slope of tob=15.69° ± 0.66°. This large difference 
in ms obtained from observations might occur because weaker active regions, 
which may not appear in white-light sunspot-group images, can be identified in 
magnetograms. However, it may also be the result of selection effects employed 
by the authors. For both case AB and case RD flux tubes, the values that we 
obtain for ?ng - j l re i ^oo -i small ii comgare^ - to - ^he i jna£netograin i data^£Stenflo and 
Kosovichev ( 20121 ), but to o large for most magnetic f ield s trengths compared to 
the sunspot-group data of lFisher. Fan, and Howard! ( 19951 ). 

The tilt-angle slopes (both mx and ttjb , as their behaviors mirror each other) 
for case RD flux tubes peak at 60 kG, whereas they peak at 40 kG for case AB 
flux tubes of 10 22 Mx. In both cases, this peak corresponds to the magnetic field 
strength where the flux tube takes the longest time to emerge (see Figure [7]). 
This behavior can be understood briefly as follows (see also Article 1). A larger 
average emergence time occurs when the effects due to magnetic buoyancy and 
the average convective downflows are of similar magnitudes, i.e. when the two 
sides of Equation 0 are approximately equal. Figure [5] illustrates this point 
well for a case AB 40 kG flux tube and a case RD 60 kG flux tube. Since these flux 
tubes spend a longer time on average in the convection zone proper, the joint 
effects of the Coriolis force and the average kinetic helicity of the convection 
simulation on the rising flux tube apex help to boost the tilt angle of the flux 
tube toward the equator, increasing the Joy’s Law trend. 

In comparison to the slopes m\ and TOb for case AB flux tubes of 10 22 Mx, 
the slopes of the best-fit lines that we obtain for case RD flux tubes are larger for 
flux tubes of 60 — 100 kG, and smaller for flux tubes of 40 kG. These differences 
are statistically significant {i.e. larger than the uncertainties). We expected that 
the overall tilt-angle best-fit slope would be reduced for case RD flux tubes 
because the value of ST = T e — T is decreased, initiating a converging parallel 
flow of mass eleinenteaDthe - ^u3^ube_a2e3^ - at^reater_de2ths_[e i- 2^Fan, Fisher, 
and McClymont, 19941 ). Due to this converging flow, the Coriolis force acts to 
tilt the flux tube apex away from the Equator in either hemisphere (opposite 
the Joy’s Law trend), thereby reducing the tilt angle. 

We suggest that the increased tilt-angle trend of 60 — 100 kG case RD flux 
tubes, as compared to case AB flux tubes, is related to a decrease in magnetic 
held strength at the flux tube apex, and that subsequently the flux tube becomes 
more susceptible to the the effects of helical convective uphows. This is capable 
of overcoming the Coriolis force induced tilt of the wrong sense {i.e. away from 
the Equator) due to a converging parallel how at the flux tube apex. However, 
this does not explain why case RD 40 kG flux tubes show a reduced tilt-angle 
trend compared to case AB hux tubes. While case RD flux tubes of <60 kG 
do not anchor, one or both of the footpoints of the buoyantly rising loops are 
usually located within the lower third of the convection zone (see Figure [4j). 
However, for case RD hux tubes of 40 kG, the hux tube as a whole resides in 
the middle of the convection zone throughout much of its evolution, where it is 
continually pummeled by convective uphows and downflows. Buoyant loops of 
these hux tubes that reach the simulation upper boundary tend to have their 
footpoints located in the middle of the convection zone (see again Figure [3J. 
Perhaps the fact that in this regime the hux tube footpoints are higher in the 
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convection zone also contributes to the reduced tilt angle trend. If stretching 
of the flux tube loop contributes to the tilting motion supplied by the Coriolis 
force ( e.g . D’Silva and Choudhuril . 1993), then loops with troughs that do not 
extend to the base of the convection zone will experience a reduced tilt angle. 
The combined effects of the reduced stretching of the loop and the increased 
converging parallel flow at the flux tube apex due to a reduction in ST may 
be strong enough at 40 kG magnetic field strengths to overcome the increased 
coupling to convective helical upflows in the non-adiabatic regime. 


4-2.2. Tilt Angle Scatter About the Joy’s Law Trend 


To quantify the scatter of the tilt angles around the best-fit line, a phenomenon 
introduced by convection, we calculate the standard deviation of the tilt about 
its fitted value following: 


E^rfo-qfit ) 2 

N 


CTfit = 


(15) 


where on is the ith tilt angle, agt is the «th tilt angle as a result of the fit 
following Method 1, N is the number of points considered, and a is in units of 
degrees. We evaluate o-fit for each field strength, where the results are given in 
columns 6 and 7 of Table [2] When calculating o-fit, we have shifted the tilt angles 
appropriately such that all tilts angles fall within the range [-90°, 90°]. This is 
unlike what was done in Article 2 where tilt angles were considered to fall within 
the range of [-180°, 180°], thereby retaining anti-Hale tilt-angle information. Due 
to the tilt-angle scatter introduced by the increased sensitivity of case RD flux 
tubes to convection, the values of Ofi t for the 40 — 60 kG case RD flux tubes 
are larger than case AB flux tubes of the sam e magnetic field strength. Using 
Mount Wilson white-light sunspot-group data, iFisher. Fan, and Howardl (119951 ) 
found that cta t ~ 30°. A value of Ofi t < 30° is found for case RD flux tubes of 
>50 kG, and >40 kG for case AB flux tubes. Similar values of erg t < 30° are 
obtained when the fit following Method 2 is used instead. 


4-2.3. Preferred Tilt Angle 

The distribution of tilt angles from our 10 22 Mx flux tube simulations are shown 
in Figure [9l where we compare both case AB and case RD flux tubes. Unlike 
the previous analyses in Section 14.21 we allow the tilt angles to fall within the 
[-180°, 180°] range: however, we only plot the range between [-25°, 55°] so as to 
highlight the difference in the two distributions near their centers. A non-linear 
least squares Gaussian fit to the distribution gives a center, or preferred tilt 
angle, of 12.7° ± 3.6° for case RD flux tubes, and 9.1° ± 3.2° for case AB flux 
tubes. The uncertainty on the preferred tilt angle is the standard deviation of the 
Gaussian fit. We note an increase in the center of the Gaussian fit for flux tubes 
with the inclusion of radiative heating, and also a wider distribution around 
the center. This increase in tilt angle reflects the results found in Section 14.2.1 1 
However, given the large standard deviation of the Gaussian fits, the centers of 
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Tilt Angle (degrees) 

Figure 9. Distribution of tilt angles in 2.5° bins, shown for both the case AB (adiabatic, 
black solid lines) and case RD (radiative diffusion, red dashed lines) simulations. For each 
simulation type, the simulation data is distinguished by diamond symbols, and the Gaussian 
fit to the data is distinguished by asterisks. The dotted line shows the zero mark along the 
horizontal axis, while the dash-dotted lines show the centers of the Gaussian fits, which are 
12.7° for the case RD simulations, and 9.1° for the case AB simulations. The addition of 
radiative diffusion to the simulation model increases the preferred tilt angle, but decreases the 
peak of the distribution. 


the di stributions are not statistically different. Using a similar analysis. iHoward 
( 19961 ') found a tilt-angle distribution that peaks between 2.5°-5° for sunspot- 
group tilt angles derived from Mount Wilson white-light photographs (1917 
1985), and between 7.5°-10° for sunspot groups derived from plage in Mount 
Wilson daily magnetograms (1967-1995). 

Although both the case AB and case RD simulations each are comprised of 
~2300 simulations, we note that the peak (i.e. height) of the distribution for 
case AB flux tubes is larger than for the RD case. There are two reasons for this 
behavior. First, the case RD distribution is broader than for case AB. Secondly, 
and most importantly, the case RD simulations have a larger percentage of tilt 
angles that are anti-Hale (> |90°|); 9.9% as compared to 5.0% for the case AB 
simulations. In comp arison, ss4 % of medium to large-sized active regions exh ibit 
anti-Hale tilt angles ( Wang and Zirinl . Il989l : Stenflo and Kosovichevl . 17(1121 ). So 
that the case RD simulations exhibit only «4 % anti-Hale tilt angles, we would 
have to exclude flux tubes with magnetic field strengths of <40 kG. In this case, 
the center of the tilt-angle distribution drops to 11.5° ± 3.3°. 


5. A Note on Flux Tube Storage in the Convective Overshoot 
Region 

In one solar dynamo paradigm, it is suggested that the dynamo generated mag¬ 
netic field is stored in a region of overshooting convective motions (overshoot 
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region) at the base of the convection zone until becoming buoyant enough to 
rise toward th e surface (e.a. ISpiegel and Weiss . 198t)l : Galloway and Weiss . 1981 


|m ieschl . ImkI i . It has been recognized that while the sub-adiabatic stratification 
of the overshoot region stabilizes magnetic flux tubes against buoyancy instabil¬ 
ities, the inflow of heat to the flux tube due to radiative diffusion forces the flux 
tube to rise quasi-statically (i.e. moving through a sequence of equilibria such 
that all forces closely balance) out of the overshoot re gion in gal year o r less , 
short compared to the wll year solar cycle period ( e.g. van Ballegooiienl . 1982 


Fan and Fisherl . 19961 : iRempell . l2003h . In this section, we review the nature of 


the solar overshoot region ISection 15. II) . and examine the upward drift of flux 
tubes in the overshoot region due to radiative diffusion (Section 15.21) . 

5.1. The Nature of the Solar Convective Overshoot Region 


Generally, the base of the convection zone is defined as the radius where the 
stratification changes from nearl y a diaba tic (V c ~ V a d) to substantially sub- 
adiabatic (V 0 < V a d) (e.y. iMieschl . 120051 ). Information about the sound-speed 
gradient in the solar interior, related to the temperature gradient, can be ob¬ 
tained from inversions of helioseismic data, revealing that the base of the con¬ 
vection zo ne occurs at 0.713 ± 0.003 I Christensen-Dalsgaard. Gough, and 
Thompson. 199it ). In reality, this observationally derived extent of the convection 
zone may include near-adiabatic portions of the overshoot region that are weakly 
sub-adiabatic (i.e. 8 = V e — V a d is negative, but V e is still very close to V a d). 
Hclioscismology has also been use d to identify the tachocline at a radius of 


0.693 ± 0.003 R 0 near the equator dCharbonneau et ali . 119991 : Basu and Antial . 
12003)1 ). indicating that the tachocline lies below the base of the convection zone, 
although it may still be located within the overshoot region. 

Convective plumes will penetrate some distance into the stably stratified 
region below the base of the convection zone before buoyancy decelerates them, 
influencing the stiffness of the transition from super-adiabatic to sub-adiabatic 
temperat ure gradient s and the thickness of the overshooting layer. As sug¬ 
gested bv iZahn (1991), efficient overshooting convection renders the temperature 
gradient weakly sub-adiabatic. Only when the convective motions have been 
decelerated to the point where turbulent mixing becomes inefficient compared 
to thermal diffusion will the temperature gradient adjust over some distance to 
that of the radiative interior. 

Much uncertainty still exists regarding the depth of the overshoot region and 
its mean sub-adiabatic temperature gradient. Helioseismic inversions find an 
upper limit on the thickness of ftdhOfjHjj. less t han 0.01 R 0 (Basu, A ntia, and 


Nara simha, 1994 Basul . 19971 : Monteiro, Christensen-Dalsgaard. and Thompson 
1994 . Early attempts to model the overshoot region predicted thicknesses of 

-1CT 6 to -10" 5 , 


0.2 H p to 0.4 Up, weakly sub-adiabatic stratifications of 8 


and a sharp transition in the temperature gradient between th e nearly adiabatic 


overshoot region and the radiative interior underneath ( e.q. 

van Ballegooiicn. 

19821: Schmitt. Rosner. and Bohnl. 1984 Pidatella and Stixl. 

19861: Skalev and 


Stix,|l99l|). However, such small values of 8 are not sufficient enough to facilitate 
the storage of magnetic flux on solar-cycle timescales (see Section E2D. It may 
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be the case that strong magnetic fields in the overshoot region can suppress 
penetrative convective motions, reducing the convective energy transport and 
convective heat conductivity, leading to a stronger sub-adiabaticity of 5 < —10 -4 
(lRemnelll2003h . 

Nu merical convection simulations with overshoot bv Brummell. Clune, and 
Toomre (l2002h exhibit a strong sub-adiabatic overshoot layer with a smoother 
transition in the temperature gradient toward the radiative value, which may be 
more hdme^ifh^l^jctu^condrtionsJrmth^Sun^ej^Christensen-Dalsgaard 
et ai, 201 il l. Th is model a lso produces larger overshooting penetration depths of 


0.4 — 2H p . Rempcll 1 2004 1 points out that the strongly sub-adiabatic overshoot 


found by numerical simulations is likely due to the larger thermal diffusivities 
required by numerical constraints, resulting in a larger energy flux. The effi¬ 
ciency of the mixing between convective upflows and downflows, as well as their 
velocities, may determine the depth of the overshoot region and the temperature 
gradient transition from nearly adiabatic to that of the radiative interior. 

For the purposes of our simulations, we have identified the base of the con¬ 
vection zone from our reference stellar structure model (Model S) as r cz b = 
0.723 Rq, the radius where the stratification changes from super-adiabatic (S > 
0) to sub-adiabatic (<5 < 0). We model a simple convective overshoot region 
where 5 decreases from 0 to —10~ 3 over a distance of ss 10 9 cm (0.01 Rq) below 
r cz b, then remains fixed at —10” 3 to the bottom of our computation domain at 
0.69 Rq. This value of <5 was chosen because it is found in Section l52?l that, S in 
the range of —10” 3 to —10 -4 is necessary for active-region-scale magnetic flux 
tubes to remain stored in the overshoot region for solar-cycle timescales. We do 
not include a transition from the sub-adiabatic stratification of the overshoot 
region to the temperature gradient of the radiative interior, as our flux tubes 
never reach such depths. 

5.2. Upward Drift of Flux Tubes in the Convective Overshoot Region 


Flux tubes initially in a neutrally buoyant state evolve very differently in the 
sub-adiabatic overshoot region compared to the super-adiabatic convection zone 
as discussed in Sections El El In the overshoot region, the buoyancy of the flux 
tube will be suppressed as the adiabatic component of the flux tube’s buoyancy 
[dAp/df] due to adiabatic expansion of the rising flux tube (Equation (1131) 1 
takes on a negative value because of the sub-adiabatic temperature gradient 
(i.e. <5 < 0). A flux tube will rise quasi-statically through the overshoot region 
if the thermal timescale over which radiative heating significantly increases the 
flux tube’s buoyancy from a neutrally buoyant state (r t h ~ 10 7 seconds) is long 
compared to the dynamic tim escale characterized by the magne tic Brunt-Vaisala 
frequency (tb ~ 10 5 seconds) ( Fan and Fisheij . ll996l ; lFanl . l2009lf . This means that 


the flux tube will evolve through a series of mechanical equilibrium states during 
_whiclMhe buoyanc y rem ains close to zero. 

Fan and Fisher ( 1996 1 derived an expression for the quasi-static rise velocity 
v T (rise speed) of the flux tube in the overshoot region, and subsequently a 
characteristic rise time for an upward drift of the flux tube out of the overshoot 
region. For simplicity, they consider a uniform horizontal flux tube that remains 
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Table 3. The rise speed v T (Column 
2) and rise time r (Column 3) for a 
100 kG, 10 22 Mx flux tube to escape 
from the overshoot region for varying 
estimates of the super-adiabaticity 
S (Column 1). In comparison, the 
length of the solar cycle is only Pill 
years. 


<5 

v T [cm s - 1 ] 

T 

—10 -3 

1.2 

26 years 

—10 -4 

12 

2.6 years 

-10“ 5 

110 

98 days 

-10-® 

670 

17 days 


nearly neutrally buoyant ( i.e. in quasi-equilibrium) in a plane-parallel overshoot 
region. Under these assumptions, the upward drift of the flux tube is found to 
be: 


v r 


Pe 


-V ad 


d Q 
d t 





(16) 


where S = V 0 — V ad , and is referred to as the sub-adiabaticity in a stably 
stratified region where 5 takes on a negative value, and 1//3 = B 2 /8Trp e . For 
typical values in the overshoot region, H p « 6xl0 9 cm, p e « 6xl0 13 g cm -1 
s -2 , (d<3/dt)i « 30 erg cm -3 s _1 , V ad « 0.4, and 7 « 5/3. Then an estimate 
for the rise speed [cm s _1 ] of the flux tube in the overshoot region becomes: 


v r 


1.2 x 10 -3 



(17) 


In the overshoot region, /3 ranges from 1.5 xlO 5 for 100 kG magnetic fields 
to 6.7xlO 6 for 15 kG magnetic fields. The rise speed of the flux tube is highly 
dependent on the super-adiabaticity [<5] o f the overshoot regi on, which is most 
likely in the range of —10 -3 to —10 -6 le.p. lRempell . 20031 200J). As <5 approaches 
zero, the term 0.12//3 becomes increasingly important. The quantity 0.12//3 is 
largest for larger magnetic field strengths; thus 100 kG flux tubes have a slower 
upward drift out of the overshoot region. Using B 0 = 100 kG, and for different 
estimates of the value 5 , the rise speed v T and time r for a flux tube to escape 
from the overshoot region are shown in Table [3] Here we have assumed the depth 
of the overshoot region to be 10 9 cm and that S is constant throughout this layer. 
For flux tubes to remain stored in the overshoot region for the length of the solar 
cycle , S needs to be somewhere between —10~ 4 and —10 ~ 3 . This is in agreement 
with iMoreno-Insertis. Schtissler. and Ferriz-Mad (1992), where it is found that 
S < -1CT 5 for B 0 < 100 kG flux tubes to be stored in the overshoot region. 
The thickness of the overshoot layer that we have assumed is slightly thicker 
than estimates derived from helioseismic inversions, but shallower than what 
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most solar models predict. A thicker overshoot region would certainly result in 
longer escape times. Rather than assuming an overshoot region with constant <5, 
a spatially varying S that increases from an average value to 0 upward from the 
base of the overshoot region would likely decrease the escape time. 

Some work suggests that convective flows may be able to help facilitate mag¬ 
netic flux storage, rather than relying on a strong sub-adiabaticity of S < — ICR 4 
to store flux tubes in the overshoot region for the length of the solar cycle. Strong 
downflow plumes in the lower convection zone could effectively pump magnetic 
flux into the overshoot region, and may also help to keep flux from escaping i 
Tobias et all 200 ll : Brummell. Clune. and Toomre , 120021 ; iBrowning et al . , 2006 


Or, perhaps an equatorward meridional flow present in the ove rshoot region may 
be sufficient to keep flu x tubes stab ly st ored f or a solar cycle (van Ball cgooiienl . 
1982 : van Ballegooiien and Choudhuri . 1988 1. However. iRempel (12003 1 found 


that such a scenario is only possible for thin flux tubes containing less than 1 % 
of a typical sunspot magnetic flux, well below the magnetic flux considered in 
this article. 

It is as of yet not clear whether magnetic flux tubes are generated at a 
tachocline interface in the overshoot region of the Sun; however, this is the 
paradigm that we adopt for this study. For the thin flux tube simulations used 
in this work and Articles 1 and 2, it is assumed that the flux tube has already 
risen out of the overshoot region. Flux tubes with footpoints that do anchor in 
the overshoot region do so because convective motions initiate buoyancy insta¬ 
bilities, which amplify undulations, causing the flux tube troughs to penetrate 
into the stably stratified plasma. In Articles 1 and 2 where radiative diffusion 
is not considered, all flux tubes that evolve with convection have footpoints 
that penetrate into the overshoot region. However, as described in Section ITT1 
when radiative heating is considered, flux tubes of <60 kG do not anchor in the 
overshoot region, while flux tubes of >80 kG do develop anchored footpoints 
that continue to drift out of the overshoot region as the flux tube evolves in 
time. 

We have performed some simulations where the flux tube is initiated in the 
stable overshoot region and allowed to rise quasi-static ally under the influence 
of radiative heating and overshooting convection. We find that the flux tube 
drifts out of the overshoot region more or less as a whole regardless of (i) the 
initial field strength of the flux tube, (ii) initial depth below r cz b, (iii) average 
sub-adiabaticity of the overshoot region , or (iv) steepness of the transition from 
S = 0 to that of the average sub-adiabaticity of the overshoot region. Once 
the toroidal flux tube has emerged as a whole from the overshoot region, the 
evolution of the flux tube proceeds in a similar way as when we initiate the flux 
tube 2.4 Mm above the base of the convection zone, as is the procedure following 
the discussion in Section [2] In order to initiate the flux tube in the overshoot 
region such that it always remains anchored except for buoyantly rising loops, 
we would have to artificially introd uce an entropy perturb ation to the flux tube 
of large amplitude, as was done bv lFan and Fisheii (1996!), rather than allowing 
convective flows to provide the perturbations self-consistently. 

Figure [TO] shows the evolution of a 5o = 60 kG, do = 12° case RD flux 
tube that is initially placed in the overshoot region at 4.9 xlO 10 cm, «10 9 cm 
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0 (degrees) 


Figure 10. Representative example of a Bq = 60 kG case RD flux tube as it rises quasi-stat- 
ically through the sub-adiabatic overshoot region and drifts upward into the super-adiabatic 
convection zone, with the base of the convection zone (r cz b) represented by the dotted line. 
The Oq = 12° flux tube is initiated below r cz b at 4.9xl0 10 cm, and the value of <5 in the 
overshoot region decreases downward from 0 at r cz b to —10 —4 over a distance of « 10 9 cm. 
The flux tube drifts upward through the overshoot region in « 2.4 years, long compared to 
the 0.2 years it takes for the flux tube to traverse the convection zone. 


below r cz b- For this particular simulation, S decreases from 0 to —10~ 4 over a 
distance of sslO 9 cm. The flux tube takes f«2.4 years to rise through the overshoot 
region, in agreement with the values quoted in Table [3] Once it emerges from 
the overshoot region the flux tube does not anchor, rising toward the surface in 
a little under two months (0.2 years). At the top of our computation domain 
(0.97 Rq), the flux tube apex emerges at a latitude of 0 em = 31.2° with a tilt 
angle of —10.1°, the appropriate sign for the northern hemisphere. All of these 
properties are consistent with the 60 kG case RD flux tubes discussed in Sections 
Eland [4] 

6. Summary and Discussion 

In this article, we modify the energy equation of our thin flux tube model to 
include the effect of flux tube heating due to radiative diffusion, unlike our 
simulations in Articles I and 2 where the flux tube is assumed to evolve adiabat- 
ically. This allows us to study the influence of radiative diffusion, in conjunction 
with solar-like convective flows, on the dynamic evolution of active-region-scale 
magnetic flux tubes. 

As a result of the inclusion of radiative diffusion to the thin flux tube energy 
equation, we find that flux tubes of <60 kG subject to convective flows are no 
longer able to anchor in the overshoot region, unlike the flux tubes discussed 
in Articles 1 and 2. Heating of the flux tubes supplied by radiative diffusion 
especially in the lower convection zone significantly increases the buoyancy of 
the flux tube earlier in its evolution. Flux tubes of <60 kG now float away from 
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the base of the convection zone before magnetic buoyancy instabilities can set 
in to anchor the troughs of the flux tube in the overshoot region. This uniformly 
increased buoyancy of the flux tube early in its evolution results in a rise time of 
<0.2 years for all flux tubes, which is significantly less than the maximum rise 
times of adiabatically evolving flux tubes by «0.5 years. 

Unlike adiabatically evolving flux tubes of 10 22 Mx, flux tubes of <60 kG 
that evolve with the addition of radiative heating exhibit a larger poleward 
deflection, especially at mid-to-high latitudes of «15°-40°. The drift of <60 kG 
flux tubes from the base of the convection zone and subsequent lack of anchoring 
facilitates a poleward slippage, which is partially mitigated by convective flows. 
When radiative heating is included in the thin flux tube energy equation, we also 
find that all magnetic-field-strength flux tubes of 10 22 Mx still exhibit a Joy’s 
Law trend in agreement with observations. However, flux tubes with magnetic 
field strengths of <40 kG exhibit too large a tilt-angle scatter and too large 
a percentage of anti-Hale tilt angles as compared to solar observations. These 
results point toward larger magnetic field strengths of >50 kG as the progenitors 
of solar active regions, inline with the adiabatically evolving flux tube simulations 
of Articles 1 and 2. 

We also discuss the problem of flux storage in the overshoot region when the 
effect of radiative diffusion on flux tube evolution is considered, suggesting that 
the time for flux tubes to escape from the overshoot region is highly sensitive to 
the value of the super-adiabaticity [J]. When we allow flux tubes to originate in 
the overshoot region, all magnetic-field-strength flux tubes rise quasi-static ally 
through the overshoot region. Once they drift out of the overshoot region as 
a whole, the evolution of the flux tube proceeds in a similar way as when we 
initiate the flux tube 2.4 Mm above the base of the convection zone. 

In our treatment of the heat input per unit volume of the flux tube [dQ v /df], 
we only include (dQ/dt)i, the contribution from the deviation of the mean 
temperature gradient in the external plasma environment from that of radiative 
equilibrium. We have here neglected (dQ/dt)^, the contribution from radiative 
diffusion across the flux tube due to the temperature difference between the flux 
tube and the external plasma environment. This is only a reasonable approxi¬ 
mation for flux tubes in the range of 15 —100 kG for magnetic flux values of 10 22 
Mx, on the order of the largest scale active regions. If (dQ/dt) 2 were included 
and we performed simulations for flux tubes of 10 20 — 10 21 Mx, we suspect that 
the average rise times for these simulations would be slightly less than those 
reported in Section HOI Following the results of Articles 1 and 2 concerning tilt 
angle trends, we expect for smaller values of magnetic flux that the scatter of 
the tilt angle about the Joy’s Law trend will increase with decreasing flux, as 
will the percentage of anti-Hale tilt angles. 

The simulations performed in this article lead us to a better and more com¬ 
prehensive understanding of the processes involved in flux emergence on the Sun 
and solar-like stars. It is likely that if flux tubes are generated in the overshoot 
region by the dynamo process, they drift upward out of the overshoot region due 
to heating of the tube plasma as a result of the deviation of the temperature 
gradient from that of radiative equilibrium. Once out of the overshoot region, 
flux tubes then drift away from the base of the convection zone as a whole due 


SOLA: rad_heat_arxiv.tex; 30 March 2015; 0:27; p. 28 



Effects of Radiative Diffusion and Convection on Rising Flux Tubes 


to an enhanced buoyancy supplied by radiative heating, only anchoring in the 
overshoot region if undular magnetic buoyancy instabilities grow sufficiently fast. 
Only when reaching the middle to upper convection zone is the dynamic evolu¬ 
tion of the flux tube well approximated by adiabatic evolution. Our simulations 
show that flux tubes do not necessarily need to have anchored footpoints in 
the overshoot region to produce emergence properties similar to those of active 
regions on the Sun. 
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